Interactions of Sodium Salicylate with β-Cyclodextrin and an Anionic Resorcin[4]arene: Mutual Diffusion Coefficients and Computational Study

The interaction between sodium salicylate (NaSal) and the two macrocycles 5,11,17,23-tetrakissulfonatomethylene-2,8,14,20-tetra(ethyl)resorcinarene (Na4EtRA) and β-cyclodextrin (β-CD) has been studied by the determination of ternary mutual diffusion coefficients, and spectroscopic and computational techniques. The results obtained by the Job method suggest that the complex formation is given in a 1:1 ratio for all systems. The mutual diffusion coefficients and the computational experiments have shown that the β-CD-NaSal system presents an inclusion process, whereas the Na4EtRA-NaSal system forms an outer-side complex. This fact is also in line with the results obtained from the computational experiments, where the calculated solvation free energy has been found to be more negative for the Na4EtRA-NaSal complex because of the partial entry of the drug inside the Na4EtRA cavity.


Introduction
β-Cyclodextrins (β-CD) are macrocycles (cyclic oligosaccharides) made up of ringshaped monosaccharide molecules, which are produced from starch by enzymatic conversion [1]. β-CD are non-toxic, edible, non-hygroscopic, chemically stable, and easily separable macrocycles. As a consequence, among the main uses that have been found for them are for the food, pharmaceutical, and drug delivery systems, and in industries such as chemical, agricultural, and environmental [1,2]. Its use in the pharmaceutical field is due to their ability to improve drug solubility and stability [3]. In addition, they have the ability to form supramolecular complexes with various chemical compounds of biological and pharmacological interest [1][2][3][4][5].
Other macrocycles of great importance in supramolecular chemistry are Resorcinarenes (RA) [6]. These compounds can be functionalized on the lower rim with hydrocarbon chains, which allows for the hydrophobicity modification of the macrocycle. Different functional groups can be added to the upper rim, such as the sulfonate group, which increases their solubility in water and allows the interaction of these compounds with solutes of pharmaceutical interest. RAs present a cavity that, under specific conditions, can be used for the inclusion of different guests by the formation of host-guest complexes (HG) [7].
Sodium salicylate (NaSal) is a compound of great interest due to its anti-inflammatory, analgesic, and antipyretic properties [8]. Additionally, it has been found that this compound induces apoptosis in cancer cells and also necrosis [8]. Due to its size and importance at pharmacological level, the study of its complexation is of great interest. Some earlier studies on the complexation of NaSal with β-CD are collected in the literature [9,10], although with results inconsistent with each other. Indeed, Junquera et al. [9], based on speed of sound and conductivity experiments, justify the presence of an NaSal-β-CD inclusion complex, according to what is usually found for most drug-CD complexes. Conversely, Deosarkar et al. [10] deduce that the external-lateral complex must be dominating over the inclusion complex. They base this conclusion mainly on quantum mechanical calculations using the DFT (Density Functional Theory) and by considering three conformers of β-CD with different patterns of hydrogen bonding between the primary hydroxyls and the salicylate (one of them with a bucket-like structure, conformer A, from which they obtain that the inclusion complex is the stable form; and the other two, conformers B and C, with a barrel-like structures, from which they deduce that the formation of the outer-side complex is the most favored). This type of discrepancy makes it convenient to further study the structure of the NaSal-β-CD complex. Likewise, the study of the interactions that may occur between NaSal and other macrocycles such as sodium 5,11,17,23-tetrakissulfonatomethylene-2,8,14,20-tetra(ethyl)resorcinarene (Na 4 EtRA) and its comparison with those that take place in the case of the NaSal-β-CD complex becomes of great interest.
Knowledge of the structure of the macrocycle-NaSal complex is of great importance. It allows for knowing if the NaSal is totally or partially encapsulated inside the macrocycle or how the NaSal is oriented inside the macrocycle cavity, offering greater protection and greater possibilities of being applied by using different delivery transit routes towards the therapeutic target, as well as greater efficiency in drug delivery and dosage [11] or, conversely, if an outer-side complex is formed, in which case the therapeutic transit routes of this macrocycle-NaSal formulation are more limited.
In this work, the diffusion and complexation of β-CD-NaSal and Na 4 EtRA-NaSal adducts is studied by using experimental (Taylor dispersion, UV-Vis, and NMR spectroscopy) techniques. In addition, both the molecular and the electronic structure of these complexes is studied by using computational techniques. The results found indicate that in both cases, host-guest type supramolecular complexes are formed between the macrocycle (Na 4 EtRA or β-CD) and NaSal, which enables the use of these results in drug delivery systems.

Results and Discussion
2.1. Ternary Diffusion Coefficients NaSal (Component 1) + Na 4 EtRA (Component 2) Table 1 shows our data of D 11 , D 12 , D 21, and D 22 . They are the mean values of between four and six replicates measured for each concentration of NaSal (1) and Na 4 EtRA (2). As it can be ascertained, the coefficients D 11 and D 22 present standard deviation values generally lower than (±0.015 × 10 −9 m 2 s −1 ), while the cross-coefficients present higher values (±0.025 × 10 −9 m 2 s −1 ). Table 1. Experimental ternary mutual diffusion coefficients, D 11 , D 12 , D 21 , and D 22 , of aqueous NaSal (C 1 ) + Na 4 EtRA (C 2 ) solutions at T = 298.15 K and P = 101.3 kPa. From this table, it can be seen that while D 11 increases, D 22 decreases with the solute fraction of NaSal (X 1 = C 1 /(C 1 + C 2 )). At the limiting situations of X 1 = 0 and X 1 = 1, the values of D 11 correspond, respectively, to the tracer diffusion coefficient of NaSal in Na 4 EtRA solutions (D 11 = 0.979 × 10 −9 m 2 s −1 ) and the binary mutual diffusion coefficient of aqueous NaSal (D 11 = 1.096 × 10 −9 m 2 s −1 ). Regarding this last value found here for D 11 , a good agreement is observed between it and the binary diffusion coefficient value already obtained in previous works (D = 1.020 ± 0.009) [12].
Similarity, D 22 values at the limiting situations of X 1 = 0 and X 1 = 1 are the binary mutual diffusion coefficient of aqueous Na 4 EtRA (D 22 = 0.711 × 10 −9 m 2 s −1 ) and the tracer diffusion coefficient of Na 4 EtRA in aqueous NaSal solutions (D 22 = 0.499 × 10 −9 m 2 s −1 ), respectively. In this case, the excellent agreement between the first limiting value and the value of the binary diffusion coefficient obtained for aqueous solutions of Na 4 ETRA at the same concentration (that is, differences around 1.5%) can also be verified [13,14].
In relation to the behavior of the cross-diffusion coefficients, D 12 and D 21 , it is observed that these are all positive, indicating that a co-current coupled flow exists. However, it is observed that D 21 is practically zero, within the uncertainty of the method (<3%). On the other hand, as the fraction of the solute NaSal increases, D 12 also increases, reaching its maximum value for X 1 →1, illustrating the existence of significant co-current coupled flows. From the values of the D 12 /D 22 ratio, it can be seen that one mole of diffusing Na 4 EtRA co-transports at most 2.4 mol of NaSal, whereas the values of the D 21 /D 11 ratio show that one mole of diffusing NaSal can co-transport up to 0.10 mol of Na 4 EtRA. At the other limit X 1 →0, D 12 is zero since it is not possible for the Na 4 EtRA concentration gradients to produce coupled flows of NaSal in solutions in which NaSal is absent.
For a better understanding of our experimental data (Table 1), a comparison was carried out between them and those predicted from Nernst equations [15] (Figures 1  and 2). These equations are very useful for the qualitative comprehension of the main characteristics of the results. In them, both the coupled diffusion mechanism and the composition dependence of the D ik coefficients are included. This theory is well described in previous works [16]. The diffusion coefficients of the ionic species Na + , Sal − , and EtRA 4− (D 0 i ) were evaluated by using Equation (3), representing F, R, T, and λ 0 i the Faraday constant, the gas constant, the temperature, and the limiting ionic conductivity.
In the present work, the limiting ionic conductivity values for Na + , Sal − , and EtRA 4− ions and of the respective diffusion coefficients are λ 0 Na+ = 50.1 × 10 −4 S m 2 s −1 [17], λ 0 Sal− = 34.5 × 10 −4 S m 2 s −1 [18] and λ 0 EtRA 4− = 229.3 × 10 −4 S m 2 s −1 [14], and     From the analysis of these figures, it is possible to verify that our data are reasonably consistent with these theoretical values. It is worth highlighting the achievement of positive values for D 12 . A possible explanation for this fact can be given based on the different mobilities of the ions involved. The values of the mobility and the diffusion coefficient of the Na + ions are higher than those of the RS 4− ions. Therefore, it must be concluded that the existence of a Na 4 EtRA concentration gradient causes the existence of an electric field that tends to slow down the Na+ ions and accelerate the RS 4− ions, since the diffusion path must always remain electrically neutral. If there is NaSal in the solution, the electric field generated by a concentration gradient of Na 4 RS will cause a coupled flow of Sal − ions in the direction of the flow of RS 4− ions. This situation collaborates to explain the reasons for the positive values obtained for the D 12 cross-coefficient. Consequently, it can be said that there are no favorable conditions for the formation of inclusion complexes between this drug and these macromolecules. Support for this fact is given by the results found from 1 H NMR spectroscopy and computational studies which suggest that there is indeed no complexation by inclusion; however, the formation of outer-side complexes is not discarded. Table 2 summarizes our data of D 11 , D 12 , D 21 , and D 22 . They are the mean values of between four and six replicates measured for each concentration of NaSal (1) and β-CD (2). As it can be ascertained, the coefficients D 11 and D 22 present standard deviation values generally lower than (±0.015 × 10 −9 m 2 s −1 ), while the cross-coefficients present higher values (±0.025 × 10 −9 m 2 s −1 ). Table 2. Experimental ternary mutual diffusion coefficients, D 11 , D 12 , D 21 , and D 22 , of aqueous NaSal(C 1 ) + β-CD (C 2 ) solutions at T = 298.15 K and P = 101.3 kPa.  From this table, it can be verified that as the solute fraction of NaSal increases (X 1 = C 1 /(C 1 + C 2 )), D 11 and D 22 increase. In the limits X 1 = 0 and X 1 = 1, D 11 represents the tracer diffusion coefficient of NaSal in 0.0100 mol dm −3 β-CD solutions, and the binary diffusion coefficient of NaSal in aqueous solutions at 0.0100 mol dm −3 , respectively. This last value is in good agreement with the value found in the literature (3.4% deviation) [19], within the uncertainty of the method.

Ternary Diffusion Coefficients NaSal
For the same composition limits, the D 22 values represent the binary mutual diffusion coefficient of aqueous β-CD (i.e., D = 0.460 ± 0.020) and the tracer diffusion coefficient of β-CD in aqueous NaSal solutions (i.e., D = 0.470 ± 0.020), respectively.
Contrarily to the previous ternary system (NaSal + Na 4 EtRA), it is observed that for the present aqueous system (NaSal + β-CD), the D 12 coefficient values are negative, reaching the maximum at X 1 →1. This means that in this case there are significant countercurrent coupled fluxes, resulting from the binding between the Sal − anion and β-CD entities. On the other hand, D 21 is practically zero for this range of concentrations, within the uncertainty of the method. From the D 12 /D 22 ratio values, it can be seen that one mole of diffusing β-CD counter-transports at most 0.50 mol of NaSal, whereas the D 21 /D 11 ratio values show that one mole of diffusing NaSal can counter-transport up to 0.04 mol of β-CD. These results show that the β-CD concentration gradients drive significant coupled flows of NaSal and, consequently, lead to favorable conditions for the formation of inclusion complexes. This effect is less pronounced when the influence of NaSal on the diffusion of β-CD is considered since the mobility of these species is expected to be quite similar.
To support these facts, the model developed by Paduano et al., whose description is well documented in the literature [20,21], can be applied in this case. Briefly, this model assumes the formation of a 1:1 supramolecular complex (Sal-β-CD), where the association constant K a is given by the equilibrium between salicylate (Sal − ) and β-cyclodextrin (β-CD) and the complex: according to where C Sal − and C β-CD are the concentrations of free (uncomplexed) Sal − and β-CD, respectively, and C Sal-β-CD is the concentration of the Sal-β-CD complex. These concentrations are correlated with each other by the following mass balance equations: Equations (6) to (9) provide the following relationships between the experimental mutual diffusion coefficients D 11 , D 12 , D 21 , and D 22 , measured for the total solute Sal − (1) + β-CD (2), and the diffusion coefficients D Sal − , D β-CD , and D Sal-β-CD which represent the diffusion coefficients of the free salicylate anion, the free β-CD and the complex, respectively, in solution. The values for these coefficients are shown in Table 3.
The diffusion coefficient value of free Sal anion (D Sal − = 0.918 × 10 −9 m 2 s −1 ) is given in reference [19], while the diffusion coefficient value of free β-cyclodextrin molecules (D β-CD = 0.436 × 10 −9 m 2 s −1 ) is estimated from D 22 at X 1 = 0. The diffusion coefficient of the Sal-β-CD complex (D Sal-β-CD = 0.400 × 10 −9 m 2 s −1 ) is estimated from the Stokes-Einstein approximation (Equations (10) and (11)), which for a given species relates its diffusion coefficient with its effective radius, r h , and therefore to the cube root of its molecular volume: k B being the Boltzmann constant, T the absolute temperature, and η the viscosity of the solvent.
From Equations 6-9 and using a value for the association constant K a equal to 80 ± 0.1 (mol −1 dm 3 ), the mutual diffusion coefficients D 11 , D 12 , D 21 , and D 22 are estimated. As it can be seen in Figures 3 and 4, said values thus estimated show good agreement with those measured experimentally (with deviations, in general, ≤3%). Furthermore, the value of K a that is selected is quite close to that found in the literature (K a = 105 ± 15 (mol −1 dm 3 )) [9].
kB being the Boltzmann constant, T the absolute temperature, and η the viscosity of the solvent. From Equations 6-9 and using a value for the association constant Ka equal to 80 ± 0.1 (mol −1 dm 3 ), the mutual diffusion coefficients D11, D12, D21, and D22 are estimated. As it can be seen in Figures 3 and 4, said values thus estimated show good agreement with those measured experimentally (with deviations, in general, ≤3%). Furthermore, the value of Ka that is selected is quite close to that found in the literature (Ka = 105 ± 15 (mol −1 dm 3 )) [9].

Complex Formation
In order to evaluate the complex formation and its stoichiometry, the Job method was applied; the results obtained for both the βCD-NaSal and Na4EtRA-NaSal systems are shown in Figures S1 and S2 in Supplementary Materials. The Job plot for both systems shows a maximum at a mole fraction of X ~0.5, indicating that the complex formation for

Complex Formation
In order to evaluate the complex formation and its stoichiometry, the Job method was applied; the results obtained for both the βCD-NaSal and Na 4 EtRA-NaSal systems are shown in Figures S1 and S2 in Supplementary Materials. The Job plot for both systems shows a maximum at a mole fraction of X~0.5, indicating that the complex formation for both systems occurs with a 1:1 stoichiometry. This fact is also in agreement with the considerations assumed for the association constant calculation. Figure 5 shows the chemical shifts (∆δ = δ free NaSal − δ complex .) of all NaSal protons as a function of the molar rate obtained from the NMR titration. For the Na 4 EtRA-NaSal complex, all values of ∆δ increase with the molar rate (Figure 5a), presumably because of the electrostatic repulsion between the free fragments; note that both Na 4 EtRA and NaSal are anionic. In addition, all the NaSal protons are shifted considerably, which can be explained either by considering a total entry of the drug into the cavity, or by the formation of an outer-side complex. We will return to this in the computational Results section.
On the other hand, ∆δ presented negative values for H 2 , H 3 , and H 4 NaSal protons (Figure 5b), which is consistent with a total entry of the drug into the β-CD cavity, because of the electronic protection that is given in the total inclusion process. A total entry of the drug into the β-CD cavity is also consistent with the fact that all NaSal protons were perturbated.

Computational Studies
Geometry optimization calculations were performed on each of the isolated carriers β-CD and Na4EtRA, with the NaSal, and the HG complexes to evaluate the complexation stability. For the HG complexes, different orientations of the NaSal guest in the cavity of β-CD and Na4EtRA were studied in such a way that both inclusion and outer-side com-

Computational Studies
Geometry optimization calculations were performed on each of the isolated carriers β-CD and Na 4 EtRA, with the NaSal, and the HG complexes to evaluate the complexation stability. For the HG complexes, different orientations of the NaSal guest in the cavity of β-CD and Na 4 EtRA were studied in such a way that both inclusion and outer-side complexes were considered. Because of the polar groups on the carriers, five different orientations of NaSal were evaluated for β-CD, whereas only three were evaluated for Na 4 EtRA. The optimized structures of these HG systems are depicted in Figure 6. Complexation energies (∆E C ) were calculated according to Equation (12) for all the orientations; the values obtained are listed in Table 4, Tables S1 and S2.  For the β-CD-NaSal system, the values of ∆EC show a higher stabilization for the inclusion complexes, likely because the amount of polar interactions formed inside the cavity of β-CD is higher compared to the outer-side complexes. The inclusion of NaSal inside the cavity of β-CD agrees with the deviation in the chemical shift observed in the NMR titration with respect to the isolated species. Among the different orientations, the hydroxyl side-on complex is the most stable (∆EC = −36.255 kcal·mol −1 ), since it is the one that  For the β-CD-NaSal system, the values of ∆E C show a higher stabilization for the inclusion complexes, likely because the amount of polar interactions formed inside the cavity of β-CD is higher compared to the outer-side complexes. The inclusion of NaSal inside the cavity of β-CD agrees with the deviation in the chemical shift observed in the NMR titration with respect to the isolated species. Among the different orientations, the hydroxyl side-on complex is the most stable (∆E C = −36.255 kcal mol −1 ), since it is the one that involves a greater amount of hydrogen bonds between NaSal and the hydroxyl groups of β-CD. On the other hand, the complexation energies for the Na 4 EtRA-NaSal system show a greater stability for the outer-side complex (∆E C = −20.362 kcal mol −1 ). In this system, both the guest and the host are charged species, therefore electron repulsions are not amenable with the inclusion of NaSal. Despite hydrogen bonds are not observed in the outer-side complex, long-range interactions allow their formation. Such outer lateral position of the NaSal, in the upper rim of Na 4 EtRA, is in line with the experimental results, thus the effect of the electric field generated by a concentration gradient of Na 4 EtRA is predominant, driving a coupled flow of Sal − ions in the same direction as the flux of EtRA 4− ions, and D 12 > 0. In summary, we found that the hydroxyl end-on and the upper outerside complexes are the dominant conformations of the β-CD-NaSal and Na 4 EtRA-NaSal systems, respectively.
We also obtained the thermodynamic properties of the complexation process (∆G C , ∆S C , ∆H C ) through analytical frequency calculations on the most stable complex of the β-CD-NaSal and Na 4 EtRA-NaSal systems (Table 4). ∆G C is negative for both systems; nevertheless, the β-CD-NaSal complex has the most negative value, which shows that the complexation of NaSal in β-CD is thermodynamically more favored. In addition, the formation of both complexes is exothermic (the values of ∆H C are −33.9 and −18.5 kcal mol −1 for the β-CD-NaSal and Na 4 EtRA-NaSal, respectively) indicating the stabilization by longrange interactions [23]. Furthermore, the higher value of ∆H C for Na 4 EtRA-NaSal can be associated with higher repulsion due to host-guest interactions.
We also calculated the solvation free energies (∆G solv ) of the complexes ( Table 4). The more negative value found for the Na 4 EtRA-NaSal complex (∆G solv = −646.8 kcal mol −1 ) is consistent with the more favorable interactions between the HG complex and the solvent. In the Na 4 EtRA-NaSal complex, the polar groups of NaSal are oriented out of the cavity of the carrier, whereas in the β-CD-NaSal system, the complete inclusion that occurs reduces the solvent-NaSal interactions. These observations can explain the negative values of D 12 for β-CD-NaSal. Most of the free salicylate is consumed by the formation of β-CD-NaSal inclusion complexes. Consequently, NaSal diffuses toward the region of higher β-CD concentration, resulting in counter-current flow to the main β-CD.
The absorption spectra of the β-CD-NaSal and Na 4 EtRA-NaSal complexes have also been simulated by TD-DFT calculations ( Figures S3 and S4, Supplementary Materials). Initially, the electronic structure of the ground state of both complexes have been analyzed in terms of their frontier molecular orbitals (Figures 7 and 8). The HOMO and LUMO of the β-CD-NaSal complex are localized on the aromatic ring of the NaSal moiety. These orbitals have a π nature, which encourages the formation of the inclusion complex (Figure 7). It is worth mentioning that in both complexes, the orbitals are strongly located on each fragment, supporting the absence of covalent interactions. Regarding the Na 4 EtRA-NaSal complex, the HOMO and LUMO are localized mainly in the RS fraction. In this complex, NaSal is stabilized by π-π stacking with the resorcinol rings of Na 4 EtRA in a lateral placement (Figure 8). Concerning the absorption wavelengths, a good agreement is obtained with the bands observed in the experimental spectra, Table 5. The band at 256 nm of β-CD-NaSal arises from the population of the A 1 state, that band is associated to an intramolecular π-π transition in the NaSal fragment. The two lowest energy bands at 210 and 182 nm result from the population of the 4 1 A and 5 1 A states, respectively. These states are also mainly the result of π-π electronic transitions in the NaSal fragment.        The bands observed in the absorption spectrum of Na 4 EtRA-NaSal at 252 and 213 are associated to the 3 1 A and 13 1 A states, respectively. These two states arise from π-π electronic transitions between the resorcinol rings of Na 4 EtRA. Nevertheless, the 3 1 A state presented a non-negligible CT transition from NaSal to Na 4 EtRA. Furthermore, the band at 213 also presented a contribution of the 16 1 A state, which is related to another CT transition between the HOMO orbital of Na 4 EtRA and the π SL+2 orbital of NaSal. The simulated UV-vis spectra of both systems, Na 4 EtRA-NaSal and β-CD-NaSal, are in line with those obtained experimentally. This fact supports the formation of the complexes described in this study.

Materials
All solutions were freshly prepared by using Millipore Milli-Q water (specific resistance = 18.2 MΩ cm, at 298.15 K). For diffusion measurements, these were prepared (at 25 • C) in calibrated volumetric flasks. For the other measurements, they were prepared by the direct weighing of both solute and solvent. To determine its concentration, the water content of the corresponding chemicals was taken into account. Before each experiment, solutions were de-aerated for 30 min, approximately.

Ternary Diffusion by Taylor's Method: Concepts and Some Experimental Aspects
For a ternary system {component (1) + component (2) + water}, the diffusion flux of component i in the direction of its concentration gradient, J i , can be determined by Fick's equations [25]: being D ii (i: 1 or 2) the main ternary mutual diffusion coefficients (representing the flow of component i caused by its own concentration gradient) and D ik the cross-diffusion coefficients (representing the coupled flow of component i caused by the gradient of the other component, k).
One of the most extensively used experimental techniques to determine values of these diffusion coefficients is the Taylor dispersion technique, which is very well described in the literature [12,[26][27][28]. In this work, each solution studied flowed in a laminar way (constant flow rate of 0.17 mL min −1 with retention times, t R ≈ 1.1 × 10 4 s, by using a Gilson Minipuls-3 metering pump) through a long capillary tube (length 3048.0 (±0.1) cm, and internal radius 0.03220 (±0.00003) cm) which was kept inside an air thermostat at 298.15 K. For given times, small amounts (0.063 mL) of a solution with a composition slightly different from that of the solution under study were successively injected into said capillary tube (with the help of a Rheodyne valve, model 5020) and its dispersion was precisely monitored, at 5 s intervals, with the help of both a Waters 2410 digital refractometer and an Agilent 34401 A digital voltmeter was placed at the end of the aforementioned capillary tube.
The voltages measured as a function of the elution time, V(t), define peaks whose shape responds to a Gaussian distribution and whose height depends on the difference in concentration between the solution under study and the injected one. The values corresponding to at least two pairs of these peaks, for each solution studied, were fitted to Equation (14): to obtain the values of the different D ik coefficients. In equation (14) D i (i: 1 or 2) represents the eigenvalues of the matrix of the ternary D ik coefficients; V 0 , V 1 and V max are the baseline voltage, the baseline slope and the peak high, respectively; and W 1 and (1 − W 1 ) are normalized pre-exponential factors.

Complex Formation
In order to study the complex formation and its stoichiometry, the Job method was performed. The experiments were carried out using a Shimadzu UV-2450 spectrophotometer. The Job method was performed by adding an aqueous host solution to an aqueous guest solution whose concentration varied from 0 to 1.0 mol fraction. The absorbance of the solutions was measured at a wavelength of λ max = 286 nm for the NaSal-Na 4 EtRA system and λ max = 296 nm for the NaSal-β-CD.
An NMR titration was performed at 298.15 K using a Varian of 400 MHz spectrometer with a 5 mm probe. Residual solvent (HOD) at pre-saturation was used to perform the spectra; 24 k data points covering a spectral width of 8 kHz were obtained, the radiofrequency excitation pulse was 45 • and the scan repetition time was 15 s to allow for complete relaxation. The resonance of the Si-(CH3)3 signal in TSP was used as internal reference, assigning it a shift value of 0 ppm. Samples were prepared in D 2 O (Eurisotop, 99.9%) as solvent, and tracer amounts of 3-(trimethylsilyl)propionic-2,2,3,3-d4 acid sodium salt (TSP), from Sigma. The NMR titration was performed by adding a carrier aqueous solution (8mM) to sodium salicylate aqueous solution (2.9mM), during the entire titration process, no precipitate was observed [29]. The results were analyzed using the HypNMR software [13].

Computational Details
The computational studies for the complex formation energies and the absorption wavelengths were performed in the framework of the Density Function Theory (DFT). In all calculations, the CAM-B3LYP functional [30] was used as in previous works [23] together with the Ahlrichs def2-SVP basis set, as included in the ORCA 4.2.1 package [31,32]. The resolution of identity approximation was taken advantage of for the calculations of Coulomb integrals with the auxiliary basis sets def2-TZVP/C [33][34][35]. The Grimme dispersion correction was also included, by the D3BJ approximation in ORCA [36] as well as the solvent effects by the Conductor-like Polarizable Continuum Model (CPCM) [37]. Water was implicitly simulated as a solvent with values of dielectric constants and refractive index of 80.4 and 1.33, respectively. The threshold for energy convergence in the self-consistent field procedure was set at 1 × 10 −8 a.u. Geometry optimization calculations were performed on the NaSal anion, the β-CD compound, the Na 4 EtRA anion, and on its corresponding HG complex. In addition, analytical frequency calculations were performed on the optimized geometries, without obtaining imaginary frequencies. The energy of the 20 low-lying singlet states was estimated by using the TD-DFT approach [38,39]. In order to simulate the absorption spectrum of the complexes, the transition wavelengths and the corresponding oscillator strengths were calculated via the electric dipole moments. A Gaussian fitting with a linewidth of 50 nm was used for plotting the absorption bands. For the depiction of geometries and orbitals, the ChemCraft visualization tools were used [22].

Conclusions
The formation of supramolecular complexes between Na 4 EtRA-NaSal and β-CD-NaSal was studied by spectroscopic and computational techniques, as well as by determining mutual diffusion coefficients. The Job plots obtained for both systems show the formation of host-guest complexes in a 1:1 ratio. For the β-CD-NaSal complex, the analysis of the diffusion coefficients and the complexation energies obtained in the computational calculations indicates that an inclusion complex occurs, while for the Na 4 EtRA-NaSal, such inclusion complex does not occur. The effect of the electric field generated by a Na 4 EtRA concentration gradient is predominant, driving a coupled flow of Sal − ions in the same direction as the flux of EtRA 4− consistent with the formation of an outer-side complex.
The solvation free energy for the Na 4 EtRA-NaSal complex has more negative values compared to those obtained for the β-CD-NaSal complex. In the Na 4 EtRA-NaSal complex, the polar groups of NaSal are oriented away the cavity of the carrier, whereas in the β-CD-NaSal system, total inclusion reduces the solvent-NaSal interactions. This is in line with the experimental results obtained for the diffusion coefficients, and also explains the behavior of the mutual diffusion coefficients.
Finally, the simulated UV-vis spectra of both systems agree with those obtained experimentally, which supports the formation and structure of the complexes studied. Funding: This research was funded by the Coimbra Chemistry Centre, which is supported by the Fundação para a Ciência e a Tecnologia (FCT), Portuguese Agency for Scientific Research, through the programs UID/QUI/UI0313/2020 and COMPETE Programme (Operational Programme for Competitiveness). It was also funded by the Faculty of Science and Vicerrectoria de Investigaciones at the Universidad de los Andes, Colombia, through the project INV-2020-99-2009; by the Faculty of Sciences and Humanities and the Comité de Ciencia y Tecnología at the Fundación Universidad de América through the project "IHU-007-2022" and finally, Minciencias by the doctoral fellowship (Doctorado Nacional-6172).